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The role of baryon-antibaryon annihilation during the hadronic stage of a relativistic heavy ion 
collision is explored by simulating the chemical evolution of a hadron gas. Beginning with a chem- 
ically equilibrated gas at an initial temperature of 165 MeV, it is assumed that the gas expands 
according to the equations of a one-dimensional Bjorken expansion. The evolution includes both 
annihilation and regeneration of baryons. From an initial time of 5 fm/c until a final time of 20 
fm/c, the number of baryons drops by approximately half for the case where there is no net baryonic 
charge. This is consistent with measurements from the ALICE Collaboration at the LHC, where the 
p/n ratio is roughly 50% lower than what comes from hydrodynamic models that ignore baryon an- 
nihilation. When the calculations are performed without the inverse process, e.g. Nir — > pp, the loss 
of baryons was found to be significantly larger. Baryon annihilation is shown to alter the extracted 
chemical breakup temperature. Assuming that annihilation cross sections are independent of the 
strangeness and isospin of the annihilating baryon and anti-baryon, the loss of strange baryons from 
annihilation is found to be similar. 
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INTRODUCTION 

High energy heavy ion collisions provide a unique opportunity to investigate the properties of nuclear matter with 
high energy density, a few GeV/fm 3 , which is several times the energy density of a proton. At these temperatures, 
> 170 MeV, matter undergoes a transition between the hadronic phase and the quark-gluon plasma phase (QGP)[TJ[2]- 
In the prevailing view, the violently interacting produced particles validate the use of thermodynamics [51-11 0| to 
describe the chemistry and hydrodynamics [TT] to describe the evolution. A typical assumption is that the hadronic 
chemistry is described by chemical freeze-out at a single temperature. This temperature is also associate with the 
QGP-hadron transition, or hadronization. Although hadrons certainly collide a few more times after hadronization, 
many models typically assume such effects to be negligible. A variety of particle ratios are reasonably well fit by 
a single temperature and a few chemical potentials to account for non-zero conserved charges such as strangeness, 
baryon number and electric charge. Additionally, fugacities have been added to account for the inability of the system 
to produce the requisite number of up, down and strange quarks for equilibrium [T2] . These pictures are all based on 
the idea that, due to the strong interactions between quarks, all partitions into various hadron species are thermally 
sampled during hadronization, and that the evolution of chemical abundances after hadronization comes only from 
decays. 

Of course, the chemical make-up does not completely freeze-out — hadrons indeed undergo chemical reactions during 
the hadronic stage. The effects of such evolution have been studied in two contexts. First, microscopic transport 
models such as URQMD (TSJ include inelastic processes. Although these models include annihilation processes, e.g., 
pp —¥ nn |141 115j . no microscopic model has yet included the inverse process. A second class of approaches involves 
solving for the chemical rates of a chemically equilibrated gas. Given that there are hundreds of resonances, such 
models typically assign fugacities to a few numbers such as effective pion number, or the net number baryons plus 
antibaryons, and then solve for the evolution of the fugacities [TBHH?] . In Reference[T!5] it was shown that reformation 
rates, e.g. nix — > pp, were important. 

The shape of pion, kaon and proton spectra can be well described by hydrodynamic models that assume chemical 
equilibrium at the hadronization temperature. However, the normalization of the proton spectra tends to significantly 
vary from experimental measurements from PHENIX at RHIC, which tend to fall ~ 40% below model predictions 
that fit the yields of pions and kaons [20] when comparing to PHENIX spectra from RHIC [21] . Similar shortcomings 
are seen at ALICE when comparing to thermal yields [H] . Using a microscopic model, a reduction of roughly this size 
was seen [TS], but those calculations ignore the inverse rates, which according to [TH] largely cancel the annihilation 
rates. The goal of this paper is to solve for the chemical rates of a wide range of hadronic species where baryon 
annihilation and reformation are consistently included. This is similar in spirit to what was done in [19] . but in more 
detail with a larger number of resonances considered. The goal will be to see whether in a consistent calculation the 
effect of annihilations might account for a ~ 40% reduction in the p/ir ratio, and to discern whether that reduction 
differs for hyperons. 

In order to concentrate on the chemistry, we employ a very simplified model of the space-time evolution. We 
consider a one-dimensional boost-invariant expansion [22j . In such a model, all fluid elements expand identically and 
their volumes increase inversely proportional to the time as measured from the beginning of the collision. We assume 
that the hadronization temperature, 165 MeV [23], is reached at a time of 5 fm/c, and that reactions cease at 20 
fm/c to be roughly consistent with expectations for 200 A GeV Au+Au collisions. In a more realistic model, the 
hadronization temperature would be reached at different times depending on the transverse position, and transverse 
expansion would accelerate the cooling at later times. The decoupling would be gradual, and the various species 
would lose thermal contact with one another and not maintain local kinetic equilibrium. Nonetheless, this picture is 
adequate to gauge the significance of baryon annihilation. 

DESCRIPTION OF CALCULATION 

We model the hadronic stage assuming that all resonances are initially populated according to chemical equilibrium 
according to a temperature To, set by default to 165 MeV. All chemical potentials are set to zero. The list of 319 
resonances extends to masses up to 2.25 GeV. Although one should have a non-zero baryonic chemical potential 
at high RHIC and LHC energies, it will be ignored for this study given that the our goal is to roughly gauge the 
importance of annihilation effects. 

After the initial time, we assume that collisions are sufficiently frequent to preserve local kinetic equilibrium, even 
though inelastic collisions are insufficient for maintaining chemical equilibrium. Each species is assigned a chemical 
potential, //$. Each fii is initially set to zero, and if chemical equilibrium were maintained each ^ would remain zero. 
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One can derive the density and entropy contributions for each species from a grand canonical ensemble, 
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where gi and fii are the degeneracy and chemical potential of hadron species i, Ei — \/p 2 + mf where m, is the 
particle mass, and V is the hadron gas volume. Once the partition function is known, the number densities, pressure, 
energy density and entropy can be found by differentiation: 
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where -K„ is the Bessel function of order n. Bose/Fcrmi effects enter through considering the sum over k. Bose effects 
are neglected for this study by considering only the k — 1 term. By including the k > 1 terms, the density of pions 
would rise at the ~ 5% level, while the effect on other species would be negligible. 

In order to maintain chemical equilibrium as the system cools, the population of heavier particles needs to decrease 
due to the fact that the thermodynamic penalty for mass is more important at low temperature. Even pions, the 
lightest mesons, need to reduce their population through annihilation and absorption to maintain chemical equilibrium 
[l6| 124] . Since baryons are heavy, their population is most sensitive to the temperature and their populations stray 
furthest from equilibrium during the rapid cooling of the hadronic stage. 

Since we are interested in simply gauging the effects of annihilation, we consider a simplified picture of the expansion, 
where densities follow a one-dimensional boost-invariant expansion [22]. This picture neglects transverse expansion 
and assumes that longitudinal acceleration is zero. Based on the assumption above, the continuity equations can be 
easily be solved in Bjorken coordinates should all chemical rates be turned off, 
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Here, To is the proper time in the local frame of the fluid cell when the system starts to expand and is set at to= 
fm/c. The evolution equations for the chemistry are modified in the presence of chemical rates, 
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where Ri is net production/annihilation rate per unit d 4 x of particles of species i. The net rates, Ri, are all zero at 
equilibrium, but as the system cools the system loses equilibrium and the rates become non-zero (usually negative) 
to push toward equilibrium. 

Once the system loses chemical equilibrium, entropy can be generated, 



dr 



[Mr)} 



T(r)' 



(5) 



Given the entropy and number densities, one can solve for the temperature and chemical potentials. Thus, Given 
the rates Ri, Eq.s (|4][5]) allow one to solve for the evolution of the number and entropy densities, or equivalently the 
chemical potentials and temperatures, as a function of r. A multidimensional Newton's method is used to numerically 
determine \Xi and T from and s. 
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One contribution to Ri is resonances decays. For the reaction A — > a\ + 0,2 + ■ ■ ■ + a n , the contribution to the 
growth rate for resonance of type A is 

1 - e ' ) (r^ Ql ... a J n A (r). (6) 

Here (r) is the decay rate of the resonance, which is the nominal width divided by the average Lorentz time dilation 
factor, 

<r) = 77r~rl / (d 3 p/E A ) m A f(p). (7) 
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Partial widths arc taken from the particle data book |25j . The contribution from the inverse reaction, ai + - ■ ■ a n — > A, 
is accounted for by the factor exp[(— fj, A + \i\ ■ ■ ■ + /i„)/T]. This forces R A to go to zero when the chemical potentials 
balance. This channel also contributes to the growth rates for a± ■ ■ ■ a n with the same magnitude but opposite sign. 
All the known hadronic decay channels from the particle data book are included in the calculation. 

The second type of chemical reaction to be considered is that from baryon annihilation, A + B — > a\ + • • • a n . The 
contribution to the growth rate of A in this case is 

Ra{t) = (l - e (<^+^+-+^-^-MB)/T) (a A+s ^ ai ... an v iel ) n A (r)n B (r). (8) 
The average (o"i> re i) for a relativistic gas is 

(o-Urel) = ( 2 ir)®n A n B J( d3 Pa/E a )(d 3 p b /E b ) f a (Pa)fb{Pb)cr(s)^(p a ■ p b ) 2 - m%m%. (9) 
We apply a simplified form for the pp annihilation cross section that is accurate for Pi a b > 100 MeV/c [26l I27j . 



a = 67 P^ 'mb. (10) 

Without particularly good justification, we further assume all baryon-antibaryon pairs have the same annihilation 
cross section as pp. Averaged over momentum, the average (av Te i) is near 5 fm-c for pp. Other inelastic processes, 
e.g., A + N — > NN, are ignored for this study as we are focusing on baryons. Such processes can change n n at the 
level on the order of 10% [HI [28]. 

Experimental data show that the most likely number of pions in pp — ► iW is ~ 5, and to simplify the calculation 
we assume that all annihilations proceed through the 5 7r production channel. For pair annihilation with nonzero 
strangeness, we employ the minimal rule for kaon production, i.c, the number of kaons in the final state equals the 
initial strangeness. For example, if the initial strangeness is unity, then final production of annihilation will have one 
kaon. We still assume final production has five particles on average. The number of pions, n,, will be 5 — tik- 

Before displaying numerical results, we summarize the assumptions in our calculation. We consider a system with 
initial charge set to zero, and at chemical equilibrium so that \Xi = for all species. The initial temperature is set to 
165 MeV, and the initial time is set to To = 5 fm/c. The calculation proceeds until freezeout at 20 fm/c. Only decays 
through the strong interaction are considered since both electromagnetic and weak decays are negligible on these time 
scales. The calculation considered 319 resonances, with masses up to 2.25 GeV/c 2 . 



RESULTS 



Figure [T] shows the evolution of the densities of pions, kaons, protons, Lambdas, Sigma, Cascade and Omega baryons. 
Results are displayed both for when annihilations are enabled and disabled. When annihilations are disabled, the 
populations of stable particles increase due to feeding from the decays of unstable particles. Decays increase the 
number of protons by roughly 50%, mainly through delta decays. The effect of baryon annihilation is to lower the 
baryon yields by ^40%, to increase pion populations by ~10% and to increase the kaon yields by ~5%. Annihilations 
reduce the yields of strange baryons by approximately the same factor as protons. For the times immediately after 
t , annihilations play no role because the inverse process exactly cancels annihilation at equilibrium. However, within 
a few fm/c, the system is no longer at chemical equilibrium and the curves with and without annihilation diverge. 

As the system cools, chemical equilibrium is lost and the chemical potentials become non-zero. Figure [2] shows the 
evolution of the temperature and chemical potentials as a function of time. The temperatures are slightly enhanced 
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FIG. 1. Hadronic densitiies, scaled by r, for the ir and K mesons, and for the p, A, S, S, and Q baryons. In the absence 
of annihilation, the proton yield increases due to the decay of resonances like the A. By adding annihilation, all baryon 
populations fall by ~ 40%. Meson populations are modestly increased by the annihilation processes. 



by annihilation, because the process pp —> 5ir creates pions with higher average energy than the thermal value. 
The chemical potentials develop more strongly for more massive particles since their populations are more sensitive 
to temperature. The lower panel demonstrates the degree to which the inverse processes cancel annihilation. At 
the initial time, chemical populations are equilibrated and the inverse processes, e.g. 5n — > pp exactly cancel the 
annihilation processes. By the end of the reaction, the inverse processes cancel ~ 55% of the annihilations. This 
underscores the importance of considering the inverse processes. 

A common use of particle ratios is to determine the chemical freeze-out temperature. Although analyses of this 
type involve a global fit, the temperature is mainly determined by ratios between the heaviest and lightest particles 
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FIG. 2. The evolution of the temperature and chemical potentials for three species as a function of time in the upper two 
panels. The lower panel shows the fraction of the annihilation rate that is uncanceled by the inverse process. Initially, the 
system is at thermal equilibrium and annihilation is perfectly canceled by the inverse process, e.g., 5tv — > pp. By the time 
breakup occurs, the inverse process cancels ~ 55% of the annihilations. 



since they are the most sensitive to the temperature. The p/ir ratio is especially important. Since such analyses, e.g. 
[B] , ignore baryon annihilation it is instructive to see the degree to which annihilation processes distort the extraction 
of a chemical freezeout temperature. Figure [3] displays the final p/ir + ratio as a function of the initial temperature 
for calculations with and without baryon annihilation. The ALICE collaboration at the LHC 29J measured p/?r + = 
0.46 ± 10%, and is illustrated with a grey band in Fig. [3| Whereas that would correspond to a chemical freeze-out 
temperature of 145 MeV if baryon annihilation were ignored, one would infer a breakup temperature near 170 MeV 
if annihilation were considered. 



CONCLUSIONS 



Figure [T] shows that annihilation processes during the hadronic stage play a significant role in the baryo-chemistry 
of the quark-gluon plasma. Calculated baryon yields that include the effects of annihilation are reduced by ~ 40%. 
This helps explain the magnitude of the p and p spectra at the LHC relative to the pion spectra without surrendering 
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FIG. 3. The final p/ir + ratio as a function of To, the initial temperature of chemical freezeout. Calculations with and 
without baryon annihilation differ substantially. If one were to use the experimental p/iv + ratio measured by ALICE to infer 
To, including annihilation would change the inferred value from ~ 145 MeV to ~ 173 MeV. 



the assumption of chemical equilibrium at the advent of the hadronic stage. Figure [T] also displays results from 
calculations where the inverse reactions, e.g. 5ir —> pp, were not included. In this case baryon yields would have 
fallen by another 40%. These results show both the importance of including annihilations, and when doing so, to 
consistently include the inverse reactions. 

The population of strange baryons is similarly affected. There are significant uncertainties in the experimental 
measurement of strange baryons, and it is premature to state whether the yields are suppressed by amounts similar 
to what is observed for protons. If it turns out that hyperon yields are not similarly suppressed, it may inspire one 
to consider whether annihilation cross sections involving strange baryons differ from those for non-strange baryons. 

Although sufficient for obtaining an understanding of the importance of annihilation chemistry during the hadronic 
stage, the model employed here is too simplistic to compare to data at even the 10% level. The assumptions are 
numerous: hadronic equilibrium at the onset of hadronization, one-dimensional boost-invariant dynamics, and local 
kinetic equilibrium that abrubtly ends at r = 20 fm/c. Further, some of the chemistry is neglected. Inelastic collisions, 
such as AN — > NN can lower the final meson populations on the order of 10%. To improve the model to the 10% 
level, one should employ a microscopic simulation of the hadronic stage. This can account for the three-dimensional 
expansion, provide a realistic freeze-out picture, and allow for the loss of local kinetic equilibrium. However, the 
current generation of such models neglects reformation processes, as well as some of the inelastic processes that affect 
the final chemical make-up measured in experiment. For instance, although baryon annihilation processes can be 
enabled in URQMD [13], the inverse processes are not. Reactions where n > 2 particles collide are more difficult to 
simulate in microscopic models. If one wishes to increase the reliabilities of these models to the ten percent level, 
such improvements are necessary. 

Figure [l] also showed the effect of annihilation on strange baryon populations. However, these results are not 
trustworthy given the lack of knowledge of annihilation cross sections for hyperons. Whereas the nucleon-antinucleon 
cross sections does not strongly depend on isospin, e.g. pp and pn annihilation cross sections are similar [30 , little 
is known about the annihilation cross sections for hyperons, either with antinucleons or antihyperons, and one must 
resort to unconstrained model calculations [3T] or to fits to particle yields in heavy ion collisions j2S] . Without direct 



experimental information about antihyperon annihilation cross sections, interpreting the yields of strange baryons 
will be problematic. In a future study, hyperon annihilation will be studied in more detail. First, it will investigated 
whether changes in hyperon annihilation cross sections significantly affect the hyperon survival probability. Secondly, 
the process where kaons scatter inelastically off nucleons to form hyperons will be investigated. In our calculations 
we see that \xk + fJ> P > /M. + Mtt- Thus, inelastic 2 — > 2 scatterings have the potential to lower nucleon yields while 
enhancing hyperon yields. In fact, such a process would successively increase the yields of hyperons of increasing 
strangeness. Although the trend is clear, the magnitude of the effect is not yet understood. 
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